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Abstract 

In this work, we present a general non-equilibrium ensemble formalism based on the subdy- 
namic equation (SKE). The constructing procedure is to use a similarity transformation between 
Gibbsian ensemble formalism and the non-equilibrium ensemble formalism. The obtained density 
distribution is a projected one that can represent essence part of (irreversible) evolution of the den- 
sity distribution, by which a generalized reduced density distribution for the quantum canonical 
ensembles is studied and applications in Cayley tree and spin network are discussed. 
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I. INTRODUCTION 



Since Gibbs synthesized a general equilibrium statistical ensemble theory, many theorists 
have attempted to generalized the Gibbsian theory to non-equilibrium phenomena domain, 
however the status of the theory of non-equilibrium phenomena can not be said as firm as well 
established] as the Gibbian ensemble theory, although great works have done by numerous 
authors^"^. The number of references along this line of research is too numerous to cite 
them all here, we just mention three significant progresses: the relevant ensembles theory 
presented-by Zubarev, Morozov and Ropke^", the Jaynes' predictive statistical mechanics 
approach^, and the generalizedi-Gibbsian ensembles theory based on the Boltzmann kinetic 
equation presented by Chan Eu^. So far the obtained non-equilibrium statistical density 
distribution formulas for the ensembles do not satisfy the original Liouville equation. Some 
researchers for that reason believe that the Liouville equation must have an extra term which 
satisfies a set of conditions assuring its irreversibility and existence of conservation laws if the 
Gibbs ensemble theory is generalized to the non-equilibrium phenomena domain based on the 
Liouville equation. But how is it possible to find this extra term which possesses universal 
irreversible characteristic to satisfy numerous requirements from a large body of models? 
This means the efforts of many school until now have not produced a universal ensemble 
theory for non-equilibrium phenomena which is comparable to the Gibbian ensemble theory 
for equilibrium phenomena. 

In this work, we present a non-equilibrium statistical ensemble formalispi based on a 
subdynamic kinetic equation (SKE)r*poted from the Brussels- Austin school^"" W and fol- 
lowed by some up-to-date works^"^". The advantage of the scheme is that SKE inter- 
twines with the original Liouville equation by a similarity transformation. If the similarity 
transformation is non-unitary, the SKE can describe the irreversible process, otherwise, it 
describes the reversible process as an equivalent equation of the original Liouville equa- 
tion. Although thererexisti several different approaches to construct SKE, that can be found 
in some publications^"^", here considering reader may be not familiar with formalism of 
subdynamics, we try to start from an update-introduction. 
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II. SUBDYNAMICS FORMALISM 



Let a quantum system S be coupled (may be strongly) to a thermal reservoir B, Hs (t), 
Hb, and H int denote the Hamiltonian of the system S, the Hamiltonian of the thermal 
reservoir B, and the interaction between S and B, respectively. The total Hamiltonian 
H (£) of the system plus the reservoir can be expressed as Hs (t) <8> Ib + Is ® Hb + H int , 
and the corresponding quantum Schrodinger equation and Liouville equation are 

i^ = H(t)f, (1) 



and 



i^ = [H(t),p], (2) 



where p — \f) (f\ is a density operator for the total system. Then one can introduce an 
orthonormal projector Pkj with Qkj = 1 — Pkj so that 

PkjH (t) Pkj = PkjHo (t) Pkj, (3) 
PkjH (t) Qkj = PkjH int Qkj. (4) 

Then the total Hamiltonian H (t) = H s (t) + H B + H in can be expressed as 

H (t) = PH (t) P + PH(t)Q + QH (t) P + QH (t) Q, (5) 

and a corresponding projected matrix is represented as 



PkjH (t) Pkj QkjH (t) Pkj 
PkjH (t) Qkj QkjH (t) Qkj 

The eigenvalue problem can be written as 

P k3 H{t)P k] Q k] H{t)P kj \ I 4> kj \f E k3 {t) \ / <f> kj 
PkjH (t) Q k j QkjH (t) Qkj I \ <f>k'j' I \ Ek'ji (t) I \ ifk'j' 



(6) 



(7) 



giving 



P kj H (t) P k j4>kj + PkjH (t) Q k j4>k'j' = E kj (t) <pkj, (8) 
QkjH (t) Pkj<pkj + QkjH (t) Qkj<Pk'j' = E k 'j' (t) 4>k>j'- (9) 
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From Eqs.Q and (jSJ), one can solve ^'j' an d <fikj respectively as 

4> Vj > = {E V ji (t) - Q kj H (t) Qkj)' 1 QkjH (t) P k j4>kj (10) 

= Ckj(f>kj) 

<j> kj = {E kj (t) - P k3 H (t) Pkj)' 1 P kj H (t) Q kj <j> wr (11) 

= Ck'j>4>k'j'- 

Substituting Eq. (|T0|) into Eq.© and Eq. ljTT!) into Eq.(j9]) respectively gives 

(P kj H (t) P kj + P kj H (t) C k , f (t)) 4> k3 = 9 feV it) 4> k3 = E k3 <t) kj , (12) 
{QkjH (t) P kj + Q kj H (t) C kj (t)) fa? = Q kj (t) (f) k 'f = E W j, (t) (j) k > f , (13) 

where introducing a creation (destruction) correlation operator (as a type of resolvent) as 

C k j (t) = Q k jC k j (t) P kj = (E k /j> (t) — Q k jH (t) Q k j) 1 Q k jH (t) P kj , (14) 
C Vj , (t) = P kj C klj , (t) Q kj = (E kj (t) - P k ,H (t) P^y 1 P kj H (t) Qk r (15) 

This shows that the ((j) k j, fik'j') are the eigenvectors of the Qyy (t) and (E k j (t) , E k /j> (t)) are 
the eigenvalues of Qk'j' (t) and H(t). This allows one to presume that the eigenvector of 
H (t) is given by f k j with the same eigenvalue E k j (t) , 

H{t) (\M) = ( W)\M V (16) 

\ l/fc'i') / \ Z k'j' 00 Ifk'j') J 
then one can find by using Eqs.liT HTBl . 

PkjH {Pkj + Qkj) fkj J / PkjHPkj PkjHQkj i i J-kjJkj 1 (17) 

QjkH (Pkj + Qkj) fkj / \ QkjHPkj QkjHQkj 

E kj (*) 
Ek>ji (t) 

Ok'ji \ I <pkj 
Qkj I \ 4>k>j> 




This gives some interesting relations as 



Pkj fkj — 4>kj — Ckj<fik'j' — CkjQkjfkj-, (18) 
Qkj fkj = 4>k'j' — Ck'j'4>kj = Ck'j'Pkjfkji (19) 



and 



fkj — Pkj4>kj + Qkj4>k'j' — (Pkj + Ck'ji) 4>kj- (20) 

Using Eq. (fi~9~]) . by introducing U k j (t) = \f k j (£)) (f k j (t)\ as an eigen-projector of H (t), one 
can construct a Schrodinger type of kinetic equation for each projected state P k jH k j (t) \<f> (t)) 



as 



i§^Pk j n kj (t)\(i>(t)) = iP kj 



|n feJ (t)]|0(t)) + n,,(t)||0(t)) 



(21) 



= P kj {[H (t) , u kj (t)} |0 (t)) + u kj (t) H (t) |0 (*)>} 

= P kj H (t) (P kj + Q kj ) U kj (t) |0 (t)> 
= P fci # (t) (P fej + C kj (t)) U kj (t) |0 (t)> , 

where, for more generality, U k j (t) can be understood as \f k j (t)) (f k j (t) in which f k j (i) G $ 
(dense subspace) and f k j (t) G $ x (generalized dual subspace of $) are defined in a Rigged 
Hilbert space, $ C 7i C $ x . This can generally provide a Schrodinger type of subdynamics 
kinetic equation (SSKE) expressed as 



with 



where 



d 

i-Qj.<f>proj CO = © (t) 0proi (t) , 
-i-Q t <t>proj (t) = © (t) 0proi CO , 



e(i)=£flyfT(t) (fly + 
b pro j and pro j (i) are defined as 

(t)) = j2 p ^k 3 (t)\^(t)), 



(pproj CO = ^ (0 CO ^i n fcj (t) 
fcjf 



(22) 
(23) 

(24) 

(25) 
(26) 



and (£) or (£) is a solution of the original Schrodinger equation in the Rigged Hilbert 



space. Furthermore, by replacing p pro j (t) = \<pproj CO) ((pproj if) 



, and using the above SSKE, 



a Liouvillian type of kinetic equation (LSKE) can be derived by 











T. oj (*)> ) ((Pproj (t) \ + \(j)proj (*)) ( (0proj (t) | 
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(27) 



© (t) \<Pproj (t)) (4>proj (t)\ - \4>proj (0) (<Voj (*) | © (*) 
[9 (t) , p proi (*)] . 
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The construction of-S 



IE or LSKE in subspace can be related to the original Schrodinger 



or Liouville equation^'^. For instance, using the relation fl20l) one have the spectral 
representation of H (t) related to 9 (t) as 



H(t) = J2Zk 1 (t)\fk J )(fk J 



z »i (*) ( p v + c *i (*)) \M iM ( p kj + Dkj (t)) 
n (t) o (t) q- 1 (t) , 



(28) 



where Djy (t) = Cjy(t), and f2 (t) = (Pj-j + Cjy (£)). The creation operator, C v = 
Q V HP U = (D u y, creates-ithe Q^-part of IT„ from the P^-part. While = H + 



Z—QvHQv 

XHiC is called collision operator^. The physical meaning of p pro j is that it represents the 
"vacuum" part of the "dynamic" part of the original density operator By, which describes the 
essence of (irreversible) evolution of the density p in its own subspace 19 '. The-isecond order 
approximation of with respect to A corresponds to the Master equation^. Moreover, 
the Boltzmann, Pauli, and Fokker-Planck equations of kinetic theory and Brownian motion 
can also be derived by using some approximation of The creation operator and 
destruction operator can also be calculated by using operator algebra to perform several 
different approaches. One of interested approaches to obtain the recurrent formulas is 



/■±oo 

C u = i drU (r) A (C v - Q v ) H (P„ + C v ) U (-r) , 
Jo 



D,. 



drU (r) A (P u + D v ) H (P„ - D U P.) U (-r) 



(29) 
(30) 



where defining U (r) = exp (—tH ). 



III. NON-EQUILIBRIUM STATISTICAL ENSEMBLES 

A marvelous remark is that the SKE seems to have the general property to approach 
various kinetic equations or Master equations, which is beyond the original Liouville equa- 
tion. As previous mentioned, the Brussels-Austin group have developed many important 
works for SKE in last two decades and have found that SKE can intertwine with the orig- 
inal Liouville equation by a similarity operator. If the similarity operator is unitary, the 
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SKE is reversible, as an equivalent representation of Liouville equation; if the similarity 
operator is not unitary, the SKE is irreversible and the corresponding evolution is not time 
symmetric. This means that the SKE can be as an appropriate kinetic equation to describe 
the irreversible process, in which the evolution operator is non-unitary on generalized func- 
tional space which is beyond the traditional Hilbert (or Liouville) space. This motivates 
one to propose using the SKE to construct a non-equilibrium statistical ensemble theory. 
The constructing procedure may be quite simple by using the "similarity transformation 
corresponding" between Gibbsian ensembles formalism based on the Liouville equation and 
the non-equilibrium ensembles formalism based on SKE: if the Hamiltonian corresponding 
to an expectation value, then the corresponding expectation of the operator should be 

Tr (Hp) = (H) — ► Tr (QPTlp) = (0) , (31) 

thus the related entropy tends to extremum, this allows one to present (by extension) a new 
canonical ensemble distribution p(9 k ) which is "vacuum" of "dynamic part" of the original 



p(Ek), as expressed by Balescu's book [19|] 



p (E k ) = Z 1 (/3, V, N) exp {-(3E k ) — > p (9 k ) = Z 1 (f3 projl V projl N proj ) exp {-(3 woj k ) 

(32) 

with the partition functions as 

Z (p, V, N) = ^ exp {-(3E k ) — ► Z ((3 proj , V projl N proj ) = ^ ex P (-PprojOk) , (33) 

k k 

(3 = (k B Ty l — > (3 proj = (ksTprajY 1 , (34) 

where 6 k is an eigenvalue of 0, (3 pro j is extended as function of position and time. In fact, 
suppose the density distribution in quantum canonical system is given by 

n(F \ / exp(-f3E kj ) 

p (E kj ) - < = , (35) 

[E fci exp(-/3£ fcj ) 

which gives the density operator p as 

^P fe> £ tJ exp(W /tjl <36) 



L—-exp(-l3H)Y,\hi) (M 



Tr exp (-8H) 

kj 

exp (-PH) 



Tr exp {—(3H) 
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Thus using the similarity transformation Q one can obtain a projected density operator p pro j 



as 



n-i n n-i ex P (-/?#) n ,on\ 

Pproj = n P n = n Trexp{ _ pH f (37) 

exp (-n^pflQ^Hfl) 



Trexpi-tt-^tttt^Hn) 
exp (-/3 proi 0) 



Trexp (-/3 pro j®)' 

This gives a precise formula of the quantum canonical ensemble for a projected density 
operator p pro j, which can be considered as generalizing the equilibrium quantum canoni- 
cal ensembles formula to the non-equilibrium quantum canonical ensembles formula in the 
sense as (1) if the similarity operator is unitary, then the new formula is just an effective 
representation of the old equilibrium quantum canonical ensembles formula because G or 
H has the same spectral structure, (2) if the similarity operator is non-unitary, then the 
new formula is an extension of the old formula and the spectrum of 6 may appear to have 
complex spectral structure that is impossible to get from the original self-adjoint operator H 
in the Hilbert space, which represents kind of non-equilibrium quantum canonical ensembles 
formula and reflects irreversibility of the system, and (3) if the similarity operator can be 
deduced by some approximations, such as Markovian/non-markovian approximations, then 
the new formula can expose some non-equilibrium characteristics, which can not be gained 
from the equilibrium quantum ensemble formulas. 

Thus it is obvious that the preceding constructed quantum formalism for density operator 
p (9k) can be extended to the classical statistical canonical ensemble by 

P (&k) = Z 1 (Pproj, Vproj, Np ro j) exp (—Pproj^k) (38) 

with 

Z (Ppro 3 , Vproj, Nproj) = J exp (-^0) dx. (39) 

In the same way, the non-equilibrium grand canonical ensembles distribution can also be 
constructed by 

P (@k) Z (Pproj, Pproj, Vproj) GXp [ Pproj^k Pproj ^proj] , (^0) 

where the partition function is given by 

Z (Pproj, Pproj, Vproj) ^ ^ J GXp [ Pproj^k Pproj Nproj] ■ (^1) 



Furthermore, the general canonical ensembles distribution may be written by 



p = Z 1 exp 



(42) 



where the thermodynamic meanings of the parameters 7&, r& can be fixed by thermodynamic 
correspondence. Again, the physical meaning of 6, f3 pr oj, Pproj and N pro j = Q^N^ are also 
the "vacuum" of "dynamic part" of the corresponding parameters, which can be functional 
of variable of the coordinate of the system; when the k system in the ensemble tends to equi- 
librium, they tend to equilibrium 0, p and A^, respectively. We want to emphasize again 
that in the book of Balescue^ the "dynamic part" means essence part of (irreversible) evo- 
lution of the density distribution, and the "vacuum" means without correlations. His work 
and Brussels-Austin school late works seem to show that the p pro j plays an important or 
influential role in the (irreversible) evok+fcion of the system by extending it to the Rigged 
Hilbert space or Rigged Liouville spaced. Using this way can one build a corresponding 
relation between equilibrium statistical ensemble formalism and non-equilibrium statistical 
ensemble formalism? The answer is confirmed because the original Hamiltonian of the sys- 
tem has corresponding relation to the collision operator by the similarity transformation. 
Thus the dynamic variables Y are usually obtained by calculated over the non-equilibrium 
statistical distribution p which is given by the proposed non-equilibrium statistical en- 
semble formulas {32} or (HDD or solution of the SKE (ETJ), (Y) = Tr (Y p (0 fc )). If the second 
order approximation of corresponds to the Master equation, the Boltzmann equation, the 
Pauli equation, or the Fokker-Planck equation, then Tr (Y p (Ok)) should deliver the expecta- 
tion of corresponding physical value in the non-equilibrium ensembles. The Eq. (1371) can be 
as starting base to get non-equilibrium statistical ensembles formulations for irreversibility, 
as demonstration of application below. 

IV. APPLICATIONS 

The pp ro j can be a generalized reduced density operator by choosing an appropriate 
projector P defined as 

«pHw.) TrBP = PpPt (43) 



Tr B exp (-/3H B 
with 



PH P= T eXP( -f"i , Tr B H„, (44) 
Treexp {-pH B ) 
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and P + Q = 1. Then the relations can be proved by assuming Hq is diagonalized and Hi nt 
is off-diagonalized: 

PHP = PH P, PH mt Q = 0, (45) 
PHQ = PH int Q, PH Q = 0. (46) 

Thus, using Eqs.([? ]) and ([? ]), one can introduce 

QIT = CPU, (47) 

TIQ = TIPD. (48) 

This gives 

p proj = \PU(f>) ^PII0| , (49) 

which is just a kind of generalization of the reduced density operator for the open system. 
This means that a generalized Markovian (or non-markovian) equation for the generalized 
reduced density operator, p pro j, in quantum canonical ensembles can be given as the formula 
( I37p . For example a generalized Markovian equation may be derived by introducing z° — 
QH Q to replace z — QHQ in the creation operator C to cancel some memory effects of C, 

Pvr ° 3 ~ TrexpHJU-0)' 

with 

e = PHP + PHQ^-y^QHP, (51) 

where z° is an eigenvalue of free Hamiltonian H . Furthermore, a Markovian equation for the 
reduced density can be obtained by the second approximation with respect to the coupling 
number A from the above equation, 

p proj = (P + D)p(P + C) (52) 

_ exp (-PprojB) 

Trexp (-ftprojQ) 
w PpP + DpP + PpC + O (A 4 ) 

w PpP + O (A 2 ) 
exp (-(3H B ) 



Tr B exp {-/3H B 

exp (-jgproj-80 

Trexp (-/Voj©') ' 



-Tr B p 
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where 

0' = PHP + PHC, (53) 

/•±oo 

C = -iX / drU (r) QHPU (-r) (54) 



x 1 

4-QH Q' 



n 



and 

P proj = PpP + 0(\*). (55) 

As an application of the above formalism to the irreversibility, let us consider a Cayley tree 
system (immersed in a nose environment) subject to a strong interaction from environment 
or an external field. The Cayley tree is a loop-free network in which there exist three classes 
of nodes, they are (1) the root nodes, which is at origin of the tree and has connectivity m, 
(2) the nodes at interface with connectivity 1, and (3) the nodes below the interface with 
connectivity m + 1. Suppose that the network start from the root of the tree with nodes 
i — 1, and link it to m new nodes i = 2, 3, • • • , m + 1, one can indicates each node with a 
subsequent number, ti indicating the time in which it arrives in the interface. At each time 
step, one can choose a node to grow, which gives rise to m new nodes. Consequently, the 
interface of the tree grows linearly in time, and the growing node is chosen at each time 
from thffi-growing number of active nodes. In order to mimic the quenched noise of the 
medium one can assign to each node of the tree an "effective" energy 6i corresponding 
to the intermediate operator by considering the strong interaction from the environment, 
and require the higher energy nodes are more likely to grow than lower energy ones. 

In fact, the total Hamiltonian of the Cayley tree system plus environment now is H = 
H + XHi, where H = Hs + Hb, H$ is Hamiltonian of the Cayley tree system, Hb is 
Hamiltonian of the environment and H\ is supposed as an interaction part of the Hamiltonian 
H by coupling to the environment with coupling number A >> 1. This kind of system 
usually is difficult to treat using the perturbative method because the series of expansion of 
the perturbation approach is related to the power of A, which is divergent. However, using 
the above proposed formula (!37l) . a generalized reduced density operator for the Cayley tree 
system, p pro j, can be written as 



.dpproj (t) 



[©, Pproj (t)} = [H + XH X C, Pproj (*)] , (56) 



dt 
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which can give 



with 



dP 

i-Qj: = POP = PH P + X 2 PH l G Q H 1 P, (57) 



G ^T^QHQ- (58) 

By taking the Born expansion Gq = 1 + XGqHi + X 2 GqHiGqHi + • • • one gets 

P(e-H )P = XPHi i _ A 1 G , Q - P. (59) 

When A >> 1, one obtain the corresponding eigenvalues of as 

n = \ {z° n + (<p n \ HxQHqQH^ 1 \<p n )) , (60) 

where <p n is the nth eigenvector of PHqP, H± is an inverse operator of Hi, and noticing 
(<p n \ PQP \f n ) = 9 n is an eigenvalue for the open system (the Cayley tree) by P tracing out 
variables of the environment. Therefore the probability p t for the active node % with energy 
9i to grow at time t can be given by 

exp + (^| HiQHoQH^ 1 {&))) 
Pi = t—^ r, (61) 

where the model depends on the parameter j3 pro j, which, can change characteristics of the 
tree. Comparing with the formula of in original case^l 

Pi = ex P tMl ; (62 ) 

Ejeint{t) ex P {-(3 z j) ' 

it can be seen that the formula of this open system still has the similar structure as the 
original one except the shift of phase (ipj \ HiQHqQH^ 1 \(fj) and /3 — > Ssjjsz. This shows that a 
node % of the tree currently possesses an an "effective" energy 8i = z®+((pi\ HiQHqQH^ 1 \<fi), 
from the original energy z®, corresponding to a random distribution p(6). This allows the 
Cayley tree network to remain, with similar characteristic as that in the original case, if 
the interaction from environment can assign to each node, such as node i, of the tree an 
"effective" energy 9^ — z® by keeping the original rule of forming Cayley tree (as a resource 
of self-orgnization) . 

Again, consider a quantum network whose nodes are composite by (elec- 
tron) spins, EdUd, {rida = d^d a , cf+ is creation operator of the spin) with in- 
teractions as connections. The environment (or control) field, Ylk o P^ko, 
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(rika = C£ a Cka, C^ a is creation operator of the fermi particle) are composed by infi- 
nite (electron) fermis. The correlation between n d is Un^n^ = ^J2 a n dcT n dcT , where U is 
correlation energy of electrons. The interaction between the network and the environment 
is A J2 k a (C£ a d a + d^Cka) , where A is coupling number. Hence, the Hamiltonian operator 
is expressed by 

H = J2 EkC+ a C k(7 + E d4da + Un d1 n dl + A ^ (C+ d a + d+C klJ ) . (63) 

When A = 0, the free Hamiltonian of the network gives two energy levels as E d and E d + U. 
This allows the nodes of the network to be possibly in three combined status: (Eq = 0), 
1 (a =| or I, E\ a = E d ), and 2 (a =| and |, 2E d + U) occupations. Hence, following wolff 
transformation^ \ we introduce three projectors to divide the total Hilbert space as three 
subspaces: 



P = (1 - n dT ) (1 - n dl ) , 

Pi = n dT (1 - n di ) + n dl (1 - n dT ) , 

P 2 = n d ^n dl , 

where P n , n = 0, 1,2 corresponds upon 0,1, 2 occupations, with 



P l = P 

1 n 1 n, 



PnPn' — 0. 



(64) 
(65) 
(66) 

(67) 
(68) 



Suppose that the total wave function \1/ is composed by vp , fa, fa in the three subspaces, 
respectively, then the Schrodinger equation = E^f can be expressed as 



^ Hqo Hqi Hq2 ^ 



#io Hu H12 
H20 H21 H22 



(fa\ 




(fa) 




= E 


^1 


[fa) 




[fa) 



where H nn i = P n HP n >, with 



(69) 



H 02 = H 2 Q = 0, 

H 01 = PqHP, = \J2 Ct (1 - n da ) d/j , 



k,a 



H 12 = P 1 HP 2 = A C+ a d a n da . 



(70) 
(71) 

(72) 



k,a 
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When the nodes are in the 1 occupation, the local extra spins appear, which allows the 
network shows type of local magnetic effect as a self organization system. Thus cancelling 
■00 and ^2 from Eq. (lp^|) and considering Eqs. (170l -(J72J), we obtain the SKE as 



9i-0i = Etipi, 

where the intermediate operator 0i is given by 

@i = Hu + H 12 (E — H22) 1 H 2 \ + H w (E — H 00 ) 1 H t 

In the second approximation with respect to A, we have 

A 2 

H12 (E — H22) H 2 i = } — (—1) C£ a C k r a >d a n da d^,n da 

' u + tid — tj d > 



01- 



k.k'.cr.cr' 



(73) 



(74) 



(75) 



H 
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(E — H 00 ) 1 H i = 2J 



X 1 



k,k',cr,a' 



E r i — Ei. 



-1) Cl a C k i a idl, (1 - n d( T>) d a (1 - nda') ■ (76) 



Considering n d = n d ^ + n d [ = 1, and S z = \ (d^d^ — d~^di), S + = d^di, S = d^d^ } the 
Eqs.f J75|) and f!76|) become to 

A 2 



Hu (E — H22) 1 H 2 \ = 2J 



k.k'.a-.a' 



U + E d -E k 



S Z (Cfct^'t ~ C H C k'i) + (77) 



H\o (E — H 00 ) 1 # 1 = 



S 2 (C^ T C fc / T - C^Cfc/J + 



(78) 



S z C^C h >^ + S C^Cfc/j + 2 ^fc^''' 

tr 

Replacing Eqs.flTT]), ([TED, and noticing iZ" u « # = Y,k,a E k C ko C kc Eq.£SD is given by 



01 = i/o + #p + 



with 



and 



U P — / j Jkk'^ ka ^k>a ► " 2^ ^fcCT^fc'^ 



k,k' ,(7 



k.k'.cr 



H ed — — J fcfc / js^ (C^Cfc/| — C^Cfc/|) + S z C^C k ^ + s C^C k > 

k,k' 

-J ^2 (Cfct^'T ~~ C^Cwi) + S z C^C k >i + S~C£fCk'. 
k,k' 



i r ' 



(79) 
(80) 

(81) 
(82) 
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where in the strong correlation condition, U — > oo, we have 



Jkk' 



x- 



J 



2 [ E k — E d U + E d — Ey 
A 2 



(83) 



\Ed — E F \ 



On the other hand, by means of the double time delay Green function (following ref. 231 ]) 

((A (t) ; B (0» = -in (t - t>) {[A (t) , B (*')]+> , (84) 
and the Furiour transformation 



((A\B)) b 



dte Mt-f) (( A (t)-,B{1/))) 



one have 



u((A\B))^{[A,B] + ) + (([A,H] \B)) b 
u((A\B))=([AM-(([A,B\ \B)) t 



All these allows one to solve the Eq. (IT9l) to obtain 



(ridi) = — arccot 



7T 



( n di) — — arccot 



7T 



E d -E F + U (n dl ) 

r 

Ed — Ep + U (n dl ) 



(85) 

(86) 
(87) 

(88) 
(89) 



where T = V 2 p^ (Ep) represents half width, in which the energy of the system distribute 
around the resonance state Ed + U (n d i) with width 2I\ and the eigenvalue of the 6 is solved 
as a complex number 

d a = E d + U (n dl ) + iT. 

When E d < E F , E d + U > E F , and \E d + U — E F \, \E d — E F \ » T, there exist a local 
solution for the Eqs. (!88|) and (!89l) : (n d i) — > 1, (ndi) ~^ 0; ^ ne network appears to have 
local magnetic vectors as a type of self-organization structure by interaction with huge nose 
environment. Its quantum statistical distribution in the (generalized) canonical ensemble 
can be given by previous formalsim as 



and 



p(M = Z~ 1 e-^ 



Z = ^ eX P (-Pprojdda) 



(90) 



(91) 
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The irreversible generalized force and the entropy for the network system can be given by 
the above canonical partition function Z formulation, i.e. 

r—$M, («) 

and 

S = k f fn Z - Pproj-^— In z] , (93) 

which shows immediately that Y and S are complex since Z including complex eigenvalue 
9 da- This means that the entropy of this irreversible system is complex! Usually, the exten- 
sion of the Hilbert space technique can be derived the complex spectrum for the self-adjoint 
operator, which demonstrates that the evolution of the intrinsic irreversible system has two 
semigroups to, represent asymmetric time evolution, which have been discussed by many 
publications' 21 ''. What is new here, through constructing the non-equilibrium ensembles for- 
mulas, we simply reveal that the complex spectrum of the G can introduce complex entropy 
and complex generalized force, which should be a characteristic of irreversible system, and 
hard to find by using the equilibrium ensemble formalism of Gibbs. 



V. CONCLUSIONS 



In conclusions, we have proposed general non-equilibrium ensembles formalism based on 
the subdynamic equation. The constructed procedure is to use a similarity transforma- 
tion between Gibbsian ensembles formalism based on the Liouville equation and the non- 
equilibrium ensemble formalism based on SKE. The obtained density distribution formula 
is a projected one that can represent essence part of (irreversible) evolution of the density 
distribution. Using this formulation, the irreversibility of the non-equilibrium system may 
emerge naturally as its entropy becomes complex, and can be exposed by calculating its 
general reduced density distribution. 
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